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Abstract: This paper presents a method of fusion of identification (attribute) information provided by 
two types of sensors: combined primary and secondary (IFF) surveillance radars and ESMs (electronic 
support measures). In the first section, the basic taxonomy of attribute identification is adopted 
in accordance with the standards of STANAG 1241 ed. 5 and STANAG 1241 ed. 6 (draft). These 
standards provide the following basic values of the attribute identifications: FRIEND; HOSTILE; 
NEUTRAL; UNKNOWN, and additional values, namely ASSUMED FRIEND and SUSPECT. The 
basis of theoretical considerations is Dezert-Smarandache theory (DSmT) of inference. This paper 
presents and uses in practice six information-fusion rules proposed by DSmT, i.e., the proportional 
conflict redistribution rules (PCR1, PCR2, PCR3, PCR4, PCR5, and PCR6), for combining identification 
information from different ESM sensors and radars. This paper demonstrates the rules of determining 
attribute information by an ESM sensor equipped with the database of radar emitters. It is proposed 
that each signal vector sent by the ESM sensor contains an extension specifying a randomized 
identification declaration (hypothesis)—a basic belief assignment (BBA). This paper also presents 
a model for determining the basic belief assignment for a combined primary and secondary radar. 
Results of the PCR rules of sensor information combining for different scenarios of a radio electronic 
situation (deterministic and Monte Carlo) are presented in the final part of this paper. They confirm 
the legitimacy of the use of Dezert-Smarandache theory in information fusion for primary radars, 
secondary radars, and ESM sensors. 


Keywords: information fusion; Dezert-Smarandache theory (DSmT) of inference; conflict redistribution 
rules; radar emitters recognition; electronic support measures (ESMs); primary and secondary radars 


1. Introduction 


Designing systems for creating a recognized air picture in the air defense system 
requires, among other factors, the development of algorithms for combining identification 
information about detected targets from various types of sensors. The basic elements of 
the air-situation-recognition system are two types of sensors: ESM sensors and combined 
primary and secondary radars (IFF: identification friend or foe). They provide identifica- 
tion information to information-processing centers that develop a recognized air picture. 
Each air-situation information-processing system should have an attribute information set, 
specifying acceptable values for the identification information of the detected targets trans- 
mitted by sensors to information-processing centers and information produced by those 
centers. This article uses a certain interpretation of attribute identification in accordance 
with the NATO STANAG 1241 standard [1,2]. It should be noted that this is one of the 
possible interpretations of the adopted definitions. It is assumed that five identification 
classes are used: three primary and two secondary ones. Sensors can transmit identification 
information in the form of a hard decision, sometimes determined as non-randomized, 
or a soft decision, sometimes determined as a randomized decision. In this paper, it is 
assumed that the sensors send identification information to the system in a randomized 
form, i.e., in the form of a basic belief assignment on the set of identification classes. This 
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assignment determines the sensor’s belief that the detected emitter belongs to separate 
identification classes. 

Another problem that should be solved by the designers of air-situation-recognition 
systems is the choice of the method of combining information from sensors. One possible 
solution is the STANAG 4162 standard proposed by NATO. It is a standard based on 
Bayesian decision functions. Another solution is to use Dempster-Shafer reasoning [3,4]. 
Works [5,6] show that such solutions have their drawbacks. The disadvantages of Dempster— 
Shafer reasoning are also confirmed in this work for situations of high conflict. These 
disadvantages can be avoided by applying Dezert-Smarandache theory (DSmT). 

DSmT in its basic version contains five rules of proportional conflict redistribution, 
namely PCR1-PCRS. Later, the authors presented a new PCR6 rule [7], which was tested in 
this paper, alongside the earlier PCR1-PCR5 rules. 

At this point, it should be noted that DSmT is not the only development of Dempster— 
Shafer theory. DSmT provides different rules for the proportional redistribution of the 
conflict mass between the resulting BBA masses in the process of information fusion. 
Examples of this development of Dempster-Shafer theory can be found in [8-10]. In [8,9], 
it was proposed to use the negation evidence (BBAs from sensors) and Deng’s [11] entropy 
to determine BBAs after information fusion. The new BBA distribution is calculated as a 
weighted sum of the original BBA, with the weights being a function of Deng’s entropy. 
The paper [10] shows the application of the negation evidence method in fault diagnosis. 

In [12], a new risk priority model based on the belief Jensen—Shannon divergence 
measure and Deng’s entropy is proposed. In the new method, Deng’s entropy and the belief 
Jensen-Shannon divergence measure are used to model the uncertainty of risk assessments 
in the “Failure mode and effects analysis” procedure and to deal with potential conflict 
information. This allows one to calculate the appropriate weighted average probabilities 
(WAPs) value. Classic Dempster’s combination rule is used to fuse data to generate 
integrated values of the risk factors. Unfortunately, the latest publications do not compare 
the DSmT method with the proposed new solutions. The examples provided there concern 
other applications than those presented in this work. 

In [13], a generalized evidential Jensen—-Shannon (GEJS) divergence to measure the 
conflict and disparity among multiple sources of evidence. This generalization was used to 
determine the weights of the information sources. Subsequently, it was used to determine 
the results of information fusion using Dempster’s combination rule. 

In [14], an extension of Dempster-Shafer theory was presented, which received the 
name complex evidence theory. It implements complex weighted discounting multisource 
information fusion. The complex evidence theory defines complex basic belief assignments 
and a complex evidential correlation coefficient. A weighted discounting multisource 
information-fusion algorithm with complex evidential correlation coefficient improves the 
performance of expert systems based on complex evidence theory. 

Another development of Dempster-Shafer theory is the generalized quantum ev- 
idence theory [15,16]. In these papers, multisource quantum information fusion was 
presented. The papers are complemented by an example of a pattern classifier from the 
motor-rotor-fault-diagnosis domain, which confirmed the efficiency of the multisource 
quantum information-fusion algorithm. 

This paper is a continuation of another work [17] and contains new research re- 
sults obtained for new scenarios of the electronic situation, new DSmT rules, and new 
information-fusion schemes. 

Below, the substantive content of individual chapters is subsequently discussed. The 
first part of the paper presents the applied interpretation of attribute identification in 
accordance with the NATO STANAG 1241 standard. It should be noted that this is one of 
the possible interpretations of the adopted definitions. It leads to the Bayesian model of the 
basic belief assignment. 

The next part of the paper presents the mathematical form of the DSmT rules PCR1, 
PCR2, PCR3, PCR4, PCR5, and PCR6 [5,18] for two sensor inputs and PCR5 and PCR6 
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for three sensor inputs, assuming the Bayesian model of the basic belief assignment of 
the hypothesis. 

The next two sections, namely Sections 4 and 5, show how to determine the basic belief 
assignment for a combined primary and secondary (IFF) radar and ESM sensors. These 
assignments are the input information in the PCR1-PCR6 information-fusion algorithms. 
Section 4 presents a method for determining the basic belief assignments (BBAs) of airborne 
targets moving in the observation space of a combined primary and secondary (IFF) radar 
sensor. This method uses the primary radar model, taken from [17]. The result of executing 
the algorithm of this method are scenarios containing reports from BBAs. Section 5 presents 
a method for determining BBAs for airborne targets that emit electromagnetic radiation 
(airborne radars and other emitters). It requires databases of reference signals of various 
airborne emitters, equipment of airborne platforms, and the nationality of the platforms. 

Each sensor report sent to the information-fusion center contains a vector of belief mass 
for all attribute identification values. The results of the proportional conflict redistribution 
sensor information, combining rules for selected deterministic and Monte Carlo scenarios, 
are presented in Sections 6 and 7 of the paper. Section 6 presents the results of research on 
the fusion of information sent only from ESM sensors. This corresponds to a situation when 
the ESM sensors operate in a system: one master station and one or two slave stations. 
Section 7 presents the results of research on the fusion of information sent from ESM sensors 
and combined primary and secondary radars. This corresponds to a situation where one 
ESM sensor (master station) and radars cooperate with the information-processing center 
(the producer of the recognized air picture). 

Conclusions are provided at the end of the paper. They confirm the legitimacy of the 
use of DSmT in information fusion for primary radars, secondary radars, and ESM sensors. 


2. Interpretation of Attribute Identification according to STANAG 1241 


The set of possible values of attribute identifications used by sensors can be adopted 
based on standardization documents of organizations that exploit these sensors [1,2,19-21]. 
This paper assumes a basic taxonomy of identification in accordance with the draft of 
STANAG 1241 ed. 6 [2]. Other similar documents may include the following standards: 
STANAG 4420 and STANAG 1241 ed. 5, which provide the following basic values of the 
attribute identifications: 
FRIEND (F); 
HOSTILE (H); 
NEUTRAL (N); 
UNKNOWN (U). 
Each of these documents contain their own definitions of the declarations. 


The following definitions of these basic values of the attribute identification are used 
in this paper (in accordance with [2]): 


e FRIEND: an allied/coalition military track, object, or entity; a track, object, or entity, 
supporting friendly forces and belonging to an allied /coalition nation or a declared or 
recognized friendly faction or group; 

e HOSTILE: a track, object, or entity whose characteristics, behavior, or origin indi- 
cate that it belongs to opposing forces, or that it poses a threat to friendly forces or 
their mission; 

e NEUTRAL: a military or civilian track, object, or entity, neither belonging to al- 
lied /coalition military forces nor to opposing military forces, whose characteristics, 
behavior, origin, or nationality indicate that it is neither supporting nor opposing 
friendly forces or their mission; 

e UNKNOWN: an evaluated track, object, or entity, which do not meet the criteria for 
any other standard identity. 


These standards bring additional values of the attribute identification: 
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e ASSUMED FRIEND (AF); 
e SUSPECT (S). 


Attention should be paid to these two last identities contained in [1], as well as their 

definitions [2]: 

e ASSUMED FRIEND: a track, object, or entity, which is assumed to be friend or neutral 
because of its characteristics, behavior, or origin; 

e SUSPECT: a track, object, or entity whose characteristics, behavior, or origin indicate 
that it potentially belongs to opposing forces or potentially poses a threat to friendly 
forces or their mission. 


The identification definitions in [1,2] can lead to different interpretations. This paper 
adopts the interpretation that is shown by the graphical form of in Figure 1. 


; AWRY 7 
cies. N — 
IWS J 
AF i, S 
Ab >, 
“Sx” U BSS 
WY YW )/ 
SS 
F H 
QY y 


Figure 1. The interpretation of STANAG 1241 using the Venn diagram. 


3. Fusion of Information from ESM Sensors and Radars in the Information-Fusion 
Center (IFC) 


3.1. Diagram of the Process of Information Fusion for Two Sensors in the 
Information-Fusion Center 

In this work, it is assumed that ESM sensors send messages asynchronously to the 
information-fusion center. These reports contain sensor decisions regarding the identifica- 
tion of objects emitting the detected signals. The set of possible identifications is as follows: 


© = {6;,=1,...,6}, (1) 


where in the following interpretation is used: 
01: FRIEND (F); 

2: HOSTILE (H); 

3: NEUTRAL (N); 

64: ASSUMED FRIEND(AF); 

65: SUSPECT(S); 

06: UNKNOWN (VU). 


According to Figure 1, the hypotheses are mutually exclusive, i.e., 


oD 


ap. (SG 
Bey ‘o fi £7. 2) 


Each sensor with the number i (i € N) sends its decisions as so-called soft decisions, 
i.e.,as BBA measure vectors (BBA: basic belief assignment). 


m; = [m;(6,),...,7;(86)]- (3) 
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A vector of generalized BBA measures for the information-fusion center should also 
be introduced as follows: 


me = [mp (01),...,mE(O6)]- (4) 


This paper adopts the Bayesian BBA model as it has been adopted as valid in the 
STANAG 4162 standard [20]. This means that Equation (5) applies in addition to (1) and (2). 


Ee mr(0;) = 08 m(0;) = 1. (5) 


In the first case that is considered, two sensors send, asynchronously in one cycle, 
one report each, containing decisions regarding the BBAs related to the target. The IFC 
system, after receiving the report from the sensor, fuses the information contained in the 
two vectors: in the current generalized BBA vector mp = [mp(61),...,mf(06)| and in the 
BBA vector m, from sensor 1 or in the BBA vector m2 from sensor 2. 

The information-fusion procedure performed in the IFC is carried out in accordance 
with the following formula: 


my; = Re(me,mj) (i = 1or 2), (6) 


wherein m, is a vector of the generalized BBA measure determined by the Rr rule based on 
the previous generalized BBA measure vector mr and the new BBA measure vector mj sent 
by the i-th sensor. The diagram of identification information fusion from the ESM sensors 
is shown in Figure 2. 


Sensor 1 IFC 
Procedure of 


ESM sensor 
and radar 
report fusion 
for two bba 


vectors 


Figure 2. The diagram of the information-fusion process in the information-fusion center (IFC) for 
two sensors. Explanations: m; a BBA measure vector of i-th sensor; mp, a generalized BBA measure 
vector that is a part of the electronic entity record in the IFC; EER, an electronic entity record in the 
IFC database. 


In the second case that is considered, two sensors send, asynchronously in one cycle, 
one report each, containing decisions regarding the BBAs related to the target. The IFC 
system waits for reports from both sensors in one cycle, using registers. Only when both 
registers are full does the IFC system perform a fusion of the information contained in 
three vectors: BBA vector mp = [mr( 01),...,me( 06)|, BBA vector m, from sensor 1, and 
BBA vector m2 from sensor 2. It should be noted that this method has a drawback: the 
information stored in the registers lose credibility. 

In this case, the information-fusion procedure performed in the IFC is carried out in 
accordance with the following formula: 


my, = Re(me,m,,mz), (7) 


wherein m, is a vector of the generalized BBA measure determined by the Rp rule based 
on the previous generalized BBA measure vector mr and the new BBA measure vectors m 
and mz sent by both sensors. The diagram of identification information fusion from the 
ESM sensors is shown in Figure 3. 
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registers IFC 
Procedure of 


ESM sensor 
Sensor 2 and radar 

report fusion 

for three bba 


Sensor 1 


vectors 


Figure 3. The diagram of the information-fusion process in the information-fusion center (IFC) for 
two sensors and electronic entity record from the IFC database. Explanations: m;, a BBA measure 
vector of i-th sensor; mr, a generalized BBA measure vector that is a part of the electronic entity 
record in the IFC; EER, an electronic entity record in the IFC database. 


Further in this paper, the combination rules of the BBA vector from the i-th sensor and 
the generalized BBA vector in the CFI are described. 


3.2. The Rules of Combination of BBA Measures Vectors 


This section presents formulas defining various combination rules for calculating basic 
belief assignments for the system shown in Figures 2 and 3. The general forms are described 
in detail in [6,18,22]. The information-fusion rules of the DSm theory are presented below 
with the following constraints: 


the properties of a set of hypotheses are described by Formulas (1) and (2); 

for the first scheme (Figure 2), the information-fusion procedure handles two informa- 
tion inputs: on one input, reports from two ESM sensors appear alternately, while on 
the second input, electronic entity records from the IFC database appear; 

e for the second scheme (Figure 3), the information-fusion procedure handles three 
information inputs: on the first input, reports from a combined primary and secondary 
surveillance radar appear; on the second input, reports from an ESM sensor appear; 
and on the third input, electronic entity records from the IFC database appear. 


3.2.1. Dempster’s Rule 

Dempster’s rule [3,4] of the BBA measure vector m; sent by the i-th sensor and the 
generalized BBA measure vector mr s in the IFC is described for each 6; € © by the 
following formula: 


LY =1,...,6 MF (K)mi(%) 


Mp (9;) =mp (9;) = 12 ae tip (O,)m;(01) 
1=1,...6 (8) 


= mer (0j)mi(8)) _ mp; (9) 
EGE waaay = Take’ 
IZk 


where in the kz; degree of conflict is defined by the formula: 


6 6 
kp = y me (0,)m; (7) = bw yy me (Ox)m; (81) , (9) 
k=1,...,6 k=11=1 
I=1,..6 Ik 
9j00,=S 
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while 
It can be noticed that 
se me (0); (01) =, (11) 
k=1,...,6 
1=1,...,6 


ee: ee 
[14.6 ity 6 
Oj Ok=D vial KFeS 
be eb 
= YE XY me(Oq)m;(9)) + LL me(O)m;(9)) = (12) 
=m = 
iAk l=k 
6.6 é 
= YY me(O)m; (01) + LY me(Oe) mj (0) = 1 
aatet| k=1 


From (12), it follows that if 


6 
y me (6); (Ox) =1 , 1e., ss me (O,)m; (07) = 0, (13) 


If 
6 
yim £(0;)m;(O,) =0 ,Le., De ME (0;) mM; (@;) = , (14) 
k=1 k=1,...,6 
1=1,...,6 


then there is no conflict. 
mp(.) is the Dempster-Shafer fusion result if, and only if, the denominator of expres- 
sion (8) is non-zero, i.e., if the degree of conflict kf; is less than 1. 


3.2.2. The Proportional Conflict Redistribution Rule PCR1 


The PCR1 rule is the simplest and the easiest version of the proportional conflict 
redistribution rule. The concept of the PCR1 rule assumes the calculation of the total 
conflicting mass (not worrying about the partial conflicting masses). The total conflicting 
mass is redistributed to all non-empty sets of hypotheses proportionally, with respect to 
their corresponding non-empty column sum of the associated mass matrix. The PCR1 rule 
is defined for every non-empty hypothesis in the following way: 


cri(0;) 
mr (8;)= mpcri (9) = 2 ; me (O)mj(%)] + a kp = 
1=1,...) 
a.r628, (15) 
cri(9;) cri(8;) 
= mp (0;)m;(6;) + ie kpj = mr;(9;) ae - ky; 
1 1 


where Cr; (0;) is the non-zero sum of the column corresponding to the hypotheses 6; in the 
mass matrix 


_— |r 
M= | (16) 
specified by the following formula: 


where: 
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e «mz; (i= 1,2) is a row vector of the basic belief assignment masses of the i-th sensor’s 
hypotheses; 

e me is a row vector of the basic belief assignments masses of the IFC system’s hypothe- 
ses; 

e =k; is the degree of mass conflict specified by the following formula: 


6 6 

kei = > me(Oq)mj(01) = )> Yo me (6,)m,; (61), (18) 
k=1,...,6 k=11=1 
1=1,...,6 l#k 
ONO) =B 


e = dr; is the sum of all non-zero column sums of all non-empty sets, as follows: 
6 6 
dr; = ) [me(0;) + mj(9;)] = )o cri()) - (19) 
j=l j=l 


In the case from this paper, dr; = 2 because 


6 6 
> me(0;) = ) | m;(6;) =1 (20) 
j=l j=l 
In addition, 
me; (0;) = me(0;)mj(0;) (21) 


3.2.3. The Proportional Conflict Redistribution Rule PCR2 


In the PCR2 rule, the total conflicting mass kr; is distributed only to the non-empty sets 
involved in the conflict (not to all non-empty sets) and taken proportionally with respect to 
their corresponding non-empty column sum. 

A non-empty set 6, € © is considered to be involved in the conflict if there exists 
another set 6; € © that is neither included in 6; nor includes a 6; such value that 6, 6; = @ 
and mr;(0,.)) > 0. The PCR2 rule is defined for every non-empty hypothesis 6; € © in 
the following way: 


m(B}) = mocral8j) = [| E_me(x)m(@)] + C08)“ kei = 
@r6) 28, (22) 
= mp(0;)m(8;) + C(0;) ED «kp = mpi(8)) + C(B}) ED - key 
where ; _ 
C(6;) = { pets ovens in the conflict, (23) 


Formula (23) can be written differently in form (25), taking into account the definition 
of the involvement in a conflict and Formula (24) [6]: 


mej(0; 0) = mp(6;) -mj(O) + me (O) - m;(8;) (24) 


1, if3Q€O,k FA]: mr (0;) -m(0,) + me(Ox) -m;(0;) >0 


C6; ji { 0, otherwise. (25) 


cr;(0;) is the non-zero sum of the column corresponding to the hypotheses in the mass 
matrix M (16) specified by the following formula: 


cri(8;) = mp(0;) + m;(9;) (26) 
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where: 

e =m; (i= 1,2) isa row vector of the basic belief assignment; masses of the i-th sensor’s hypotheses; 
e me isarow vector of the basic belief assignment masses of the IFC system’s hypotheses; 
e = ky; is the degree of mass conflict specified by Formula (18); 

e = er; is the sum of all non-zero column sums of all non-empty sets involved in the 


conflict, as follows: 


6 6 
eri = )) [mr(8;) + mi(8))] = DL, cri(O;) = Y) C(%) [me (8;) + mi(8;)] =), C(8)) - cri(9)) (27) 
j€CF jeCE jal j=l 
where 
CF = {j =1,...,,6 : VO,E@O me; (0; Oz) > o} (28) 


and mr;(0; 9 0) is defined by (24). 
In addition 
mpj(0;) = mp(0;)m;(8;) (29) 
It is shown below that in the case of data used in numerical experiments (Section 6), 


er; = 2, which means that the PCR2 rule is equivalent to the PCR1 rule. The BBA vectors 
used there contain values less than 1, which means the following: 


Vj=1,...,6 : mp(6;) < 1A m(6;) < 1 (30) 


It follows that each BBA vector contains at least two non-zero components, that is 
dj=1,...,6, dk=1,...,6 with k F j, such that 


0< me (6;) <1A0< me (Ox) <1 (31) 


0< mj (6;) <1A0< mj (A) <1 (32) 


From (31) and (32), it follows that if me (0;) > 0, then there exists at least one value 
k # j, such that m;(60,) > 0, which can be written in the following form: 


Wi = 14a, 6 3 0 < me (6j) <1 => dk...... ,6,k Aj: 0< m(&)< 1 (33) 


From (33), it follows that 


Vj=1,...,6 : 0<mp(0) <1 > 3k=1,...,6,k#j : me(;)-mj(6) > 0 (A) 


The same applies to the following: 


Vj=1,.-.,6: 0<m,(0))< 1 > Sk=1,...,6k Aj: m;(0j))-mp(%)> 0 (35) 


J 


Taking into account (34), (35) and (25) can obtain the following: 


fea 0 < mr(6;) <1 =} ak=1,...,6,k Aj such that 


36 
me (0;) Z m; (0x) a m;(0;) : me (Ox) >0 oe 
Vi=1.%470 7 0 m;(8)) <1 = ak=1)..4,6,6 $7 such that (37) 
mj(0;) - mp(Ox) + me (Oj) - mji(Oq) > 0 
From (36) and (37), it follows that 
Vj =1,...,6: 0 < me(0;) <1 => C(@)=1 (38) 


Vj=1,...,6 : 0<m;(6)) <1 > C(a)=1 (39) 
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This means that any hypothesis with a non-zero BBA value for any of the two sensors 
is involved in a conflict. 
From (27), it follows that 


6 6 6 
eri = ) | C(8;) [me(8;) + m;(9;) | = )) C(Oj)me (Oj) + Y) C(8;)m:(9;) (40) 
j=1 j=l j=l 
Using (36), (37) and (40), the value ef; is determined. 
Because 
6 6 6 6 
Le CO)me(6j)) = DL C(0j)mr(O;) + DL) C(Oj)me(8;))= J) me(8;) =1 (41) 
j=l j=l j=l j=l 
me (0;)>0 me (0;)= mp (8;)>0 
6 6 6 6 
y= C(6;)m;(4;) = ; C(0;)m;(0;) + a C(6;)m;(9;) = yy m;(0;) =1 (42) 
=1 = j=1 =1 
: mj(0;)>0 a Se, m;(6,)>0 
we obtain 
6 
er; = )_ C(0;)mp(8;) + ) | C(0;)m;(0;) = 2 (43) 
j=l j=l 


Considering (43), it can be said that in this case, the PCR2 rule is equivalent to the 
PCR1 rule. For this reason, the results of the PCR2 rule are not presented in Section 6, as 
they would be identical to the results of the PCR1 rule as only the Bayesian BBAs are used 
in this application. 


3.2.4. The Proportional Conflict Redistribution Rule PCR3 


In the PCR3 rule, the partial conflicting masses are distributed instead of the total 
conflicting mass, kr;, to the non-empty sets involved in the partial conflict. If an intersection 
is empty, for instance 6,6; = ©, then the mass m(6,. 6;) of the partial conflict is 
transferred to the non-empty sets 6; and @ proportionally, with respect to the non-zero sum 
of masses assigned to 6; and, respectively, to 6; by BBAs mp(.) and m;(.). The PCR3 rule 
works if at least one set between 6; and 6; is non-empty and its column sum is non-zero. 

The PCR3 rule is defined for every non-empty hypothesis 6; € © in the following way: 


my (6;) = mpcr3(0;) = _ : me (Ox)m;(8) | + | cri(;) > SECRS (8,64) | = 
I=1).5.6 0,10;—2 
8,78; —8; 
me;(0;) aE: [cri(8j) (Sr (8j-9)| 
kA 
(44) 
where 
me (Ox) m; (0; ) +e (8; )m (Ox) 1D. 
SEER, 6x) _ Cri(0j) Feri (Ox) for cr; (6;) + cri (O) #0 (45) 
O for cp;{(9;) + cri(%) = 0 
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SEF (0;, Oy.) = 


cr;(6;) is the non-zero sum of the column corresponding to the hypotheses 6; in the 
mass matrix M (16), specified by the following formula: 


cri(8;) = mp(0;) + m;(9;) (46) 


3.2.5. The Proportional Conflict Redistribution Rule PCR4 


The PCR4 rule redistributes the partial conflicting masses only to the sets involved 
in the partial conflict in proportion to the non-zero mass sum assigned to 6, and 6; by the 
conjunction rule according to the following formula: 


me (8;) = mpcra(6j) = mei(0j) + mFi(8)) LE SEER), 0) = 


=e 
NO; =D 47 
= me;(0;) + mpi(6j) or OO 0) ( ) 
kA 
where 
(8,78 
saa for cri(8)) + cri(8x) #0 and mp;(6;) -mri(6,) #0 
(8:08 
aby testa for cpj(0;) + cri(O) AO and mp; (Oj) - mpj(Ox) = 0 ) 
0 for cE (0;) + cri (A) =0 
wherein 
me; (0; NO) = me (Ox) mj (0; ) + me( 0;)m; (Ox) (49) 
me;(0;) = me(0;) - m;(9;) (50) 
cri(0;) = mp(0;) + m;(0;) (51) 


If at least one of the BBAs, me(.) or m;(.), is zero, the fraction is discarded and the 
mass mf;(6; 1 0;) is transferred to 6; and 6; proportionally, with respect to their non-zero 
column sum of masses cf;(6j). 


3.2.6. The Proportional Conflict Redistribution Rule PCR5 for Two BBAs (Two Sources) 


Similar to the PCR2-PCR4 rules, PCR5 redistributes the partial conflicting mass to 
the hypothesis involved in the partial conflict. PCR5 provides the most mathematically 
precise [6,18,22] redistribution of conflicting mass to non-empty sets in accordance with the 
logic of the conjunctive rule. However, it is more difficult to implement. The PCR5 rule is 
defined for every non-empty hypothesis 6; € © in the following way: 


m',(6;) = mpcrs(6j) = mpi(8)) + Yo SEER (6,0) = mpi(6;)) + Yo Ser® (6), 6x) 


k=1,...,6 k=1,...,6 
Oj = kFj 
(52) 
where 
me (8;)?-m;(6 i __-mj(0;)?-me (8) 
PCRS = ME (6; ) +m, (6 m; (6; ) +i (8) 
Sri (89k) = 9 map(8;) + m(6x) 92 Oand (6 ;) + mp (0) #0 o) 
0 for inte (8; ) + m;(0,) = 0 or m;(8;) +me(6,) =0 
wherein 
mr; (6;) = me(8;) - m;(8;) (54) 


In Formula (52), the component SECRS is equal to zero if both denominators are equal 
to zero. In Formula (53), if a denominator is zero, then the component is discarded. 
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SPo” (Bj, 9K, 91 


Sie (6;,0,) = 


SF (8; 9) = 


3.2.7. The Proportional Conflict Redistribution Rules PCR5 and PCR6 for Three BBAs 
(Three Sources) 

In [6,22], improved proportional conflict redistribution rules of the combination of ba- 
sic belief assignments PCR6, PCR5+, and PCR6+ are presented. The authors point out that 
these rules should be applied if, and only if, we are to combine more than two BBAs. If we 
only have two BBAs to combine (s = 2), we always obtain mpcrs = Mpcrs5+ = Mpcr6 = MPCR6+, 
because in this case, the PCR5, PCR5+, PCR6, and PCR6+ rules coincide. Below are the 
formulas that define the PCR5 and PCR6 rules for three BBAs. 

The PCR5 rule for three BBAs (three sources) is defined for every non-empty hypothe- 
sis in the following way: 


m" (6;) 
m,(0)) = mpcrs(0)) = ——“—~ (55) 
ewe "YS m' (6) 
wherein 
(8) + ESE (Gj- rH) +E SUPE (Gj %)++ EL S2EF (6), Ok) = 
[=1,...6 6,0; and 
OK8)6; =O 
Lu sige (0), 84,61) + DL SUB? (8, O)++ LC S2eG°(Gj,0) = 
=1,..6 1=1,. ae PT ace (56) 
ipikizk i#k TP 
Sia (OFebh) SS1ES 0,0) + Ste be) 
YT iKiet 
__ mp(8;)? = 14 (Bq) « m2(8,) me (81) - 121(8))* + 12(OK) ME (OK) = 14 (81) - 285)? 67) 
me (0;) + my (0,) + m2(6;) — me(6)) +111 (8j) + m2(Oq) ° me (O% Ee ACEC) 
me (0;)? - my (Oq) m2 (O) — me(Oq) «114 (8;)? «my (Ox) te (Bx) » 14 (Bq) - 112 (8;) (58) 
mr (0;) + m1(0,) + m2(O,) me(O,) + my (6;) + m2(6,) me(O,) + m4(0,) + (8) 
me (0;)? - 114(0;)* = m2(Oq) mp (O) - 114 (8;)? - m2 (0;)? me (O;)? - my (8) - m2(8;)* 
F\j AY 2WVk) | MN F\UK WWW DAU f 9 MEENUG 1\Wk 2\9; (59) 
F(0j) + m1 (0j) + m2 (Oq) © me (Ox) +1 (8;) +m2(8;) © me(B;) + 111 (Bq) + m2(6;) 
m!' (0; ) 
m, (0;) = mpcrs(0)) = ——“—~ (60) 
BN] J es ,m" (6 FB) 


In Formulas (57)-(59), if a denominator is zero, then the component is discarded. 
The quotient in Formula (55) ensures the normalization of the BBA vector Mp, which 
ensures the following: 
6 
yo me (6;) = ) mpcrs(6;) = 1 
i=1 i=1 
The PCR6 rule for three BBAs (three sources) is defined for every non-empty hypothe- 
sis 0; € © in the following way: 


m" (0;) 


—— 61 
Fe, m" (8) oe 


mp (0;) = mpcrs(;) = 
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wherein 
m" (6j) = mp2 (0;) + —, SECR (4 ;, Ox, 81) + , S1ESRE (0, 0) + 
[=1,...6 8,16; 
88,8 =O 

+E SUPER (6; 4) = 

0,76;=2 

= me12(8;) ae. ty L ate Cre Ok, 61) seo oo Eras (Oy 0.) + 

k=1,...,6 1=1,..6 k=1,...,6 (62) 
kei IAIN FR i#k 
a S2P 3° (8; 8%) = 
j-k 


=men(G)+ E | 2d SECR6(8;,.0., 01) + SIPCR(0;, 6.) + S2PER6 (A, ,) 


k#j_ LIAjAZk 
= (6;)2-m, (B,)-m3(64) (6,)-m, (6;)?-ma(6) 
PCR5 : ie Mp 9; “My OK “m2 0; f Me 0) “my i “M2 OK j 
Sra” (8j- 91) mp (0;) +11 (O,)+m9(8)) mp (6) +m (6)) +19 (6) (63) 
__me (Ox): (81) -112(8;)" 
"ate (Ox) +1 (8;) +119 (8;) 
PCR5(p. __me(6;)? +1 (8,)-m2(8,) me (8-711 (8; )* 112 (8) 
Sli (6;, 9x) me (8; ) +m (8) +1112 (O) | me (8) +m (6; ) +712 (0) (64) 
tp (8x): (8x) (8))? 
: me (Of) +my (Ox) +m (8;) 
pene me (0;)? « my (0;) + my (O,) + me(O;) - m1 (0;)? - m2 (O) 
S2r12 (9), 0) = 
me (8;) + my (0;) + m2 (Ox) 
mp (8) + my (;)* » my (6;) + mp (8x) «11 (8;) - 2(8;)? (65) 
me (Ox) +My, (0;) + Mp2 (0;) 
mp (6))> - 114 (8) + mz (8;) + mp (8) - 114 (84) + 12 (8))? 
mr (0;) =e m1(Ox) a m(6;) 
and 
mpry2(0;) = mp(0;) - m1 (0;) - m2(8;) (66) 


In Formulas (63)-(65), if a denominator is zero, then the component is discarded. 
The quotient in Formula (61) ensures the normalization of the BBA vector mp, which 
ensures the following: 


Me 


6 
mr(0;) = )- mpcrs(6;) = 1 
= 


i=1 


Comparing the two fusion schemes (Figures 2 and 3), it should be noted that sequential 
and global information fusion generally produces different results [18], i-e., 


PCR5(me, m1, m2) #4 PCR5(PCR5(me, m1), m2) 4 PCR5(PCR5(me,mz),m,) (67) 


In addition, the article experimentally verified the theorem on the inequality of the 
results of both PCR5 and PCR6 rules for three BBAs (three sources) presented in [18]: 


PCR5(me, m1, m2) #4 PCR6(me, m1, m2) (68) 


4. Basic Belief Assignment for Combined Primary and Secondary Surveillance Radars 
Combined primary and secondary (IFF) radars are the main source of identification 

information regarding air and maritime objects. A primary radar only yields the detection 

of an object in a supervised area. The detection of the object is the precondition for sending 
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a request to the object by the secondary radar (interrogator). Interpretation of the object 
response is dependent on the type of request. The so-called civilian modes only yield a 
determination of whether the detected object replies to an interrogation or not. 

This paper assumes that the analyzed radar sensor consists of two radars: primary and 
secondary. Therefore, the probability of the correct detection and the correct identification 
of a target is expressed by the following formula: 


Pyj = Pa Pier (69) 


where P is the probability of correct detection of the target by a primary radar, and Pyrr is 
the probability of a correct reply to an interrogation. If a target is detected by the primary 
radar and there is a lack of proper identification by the secondary radar, it can be assumed 
that the target has a value of attribute identification of UNKNOWN—U. Thus, the following 
relation can be written: 

m(U) = Pa(1— Pree), (70) 


where m(U) is the mass of probability for a value of UNKNOWN identification attribute. 
A method for calculating the probabilities Py and Pirr is presented in [7,17,23]. 
This section explains the way the remaining mass of probability is calculated (1 — m(U)). 
It is assumed there that every simulated target should have a base value of attribute identi- 
fication from the set as follows: 


Zpr = {Np, Fp,Hp} (71) 


where: 


e Nz: base NEUTRAL identity; 
e Fg: base FRIEND identity; 
e Hz: base HOSTILE identity. 


STANAG 1241 introduces, in addition to the basic set of attribute identification val- 
ues, secondary (additional) attribute identification values: SUSPECT (S) and ASSUMED 
FRIEND (AF). According to Figure 1, a table of possible attribute value transitions between 
set (10) and the set of secondary attribute identification values can be introduced: 


ZsI — {Ns, Fs, Hs, AF, S} (72) 


The belief mass values contained in Table 1 determine how the mass of the base belief 
assignment is transformed into the mass of the secondary belief assignment. They can be 
estimated as empirical frequencies based on recorded archive events. 


Table 1. Transformation of the base belief assignment mass into the secondary belief assignment mass. 


Base Identification > Fg Ng Hp 
Fs m(Fs|Fp) 0 0 
Ns 0 m(Ngs|Np) 0 
Hs 0 0 m(Hs|Hp) 
AF m(AEF|Fg) m(AE|Np) 0 
5 0 m(S|Fz) m(S|Hp) 


Of course, the normalization conditions are satisfied: )7.<7,,™(x|Fg) = 1, Uxez,,m(x|Np) = 1, 
and )ivezo,m(x|Hp) = 1. 

The final values of the belief mass of secondary attribute identification values are 
calculated according to the formulas as follows: 
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1. For a target with the FRIEND base value of an attribute identification, 


) 
m(AE) = m(AE|Eg) (1 — m(U)); 
m(Fs) = m(Fs|Fp)(1 — m(U)). 


2. For a target with the NEUTRAL base value of an attribute identification, 


(U) = P3(1 — Pree); 
eee = m(AF|Ng)(1—m(U)); 
( 


3. For a target with the HOSTILE base value of an attribute identification, 


m(U) = Py(1 — Pier); 
m(Hs) = m(Hs|Hg)(1—m(U)); 
m(S) = m(S|Hg)(1— m(U)). 


Other final values of the belief mass of secondary attribute identification values are 
equal to zero. 


5. Basic Belief Assignment for ESM Sensors 


ESM sensors consist of passive receivers and direction finders, which allow them to 
capture emitter signals coming from certain directions. In this way, the electronic recog- 
nition system can receive, among other data, information on radar emitters mounted 
on air or maritime platforms. Reports sent from the ESM sensors include, among oth- 
ers, the characteristics of the intercepted signal, the emitter’s azimuth, and the so-called 
identification information. 

This paper also assumes that sensors are equipped with specialized databases called 
the databases of emitter signal patterns, in which information about previously captured, 
processed, analyzed, recognized, and described radar emitter signals is stored, along with 
additional information about the type and mode of the emitter work, the platform on 
which these emitters can be installed, and the national or organizational affiliation of these 
platforms. The detected signals are the subject of an analysis procedure, which yields 
the determination of the so-called distinctive features of the signal, and then assigns this 
information to a specific electronic entity (already existing or created ad hoc) [24]. The 
basis for assigning distinctive information to an electronic entity is the azimuth angle of the 
incoming signal. 

In the case of a high density of targets, identification information may fluctuate due to 
incorrect assignment of signal information to the electronic entity [25]. The impact of this 
negative phenomenon can be significantly reduced by an efficient estimation of the emitter 
positions [24]. Assuming that the sensors send all reports on the tracked electronic entities 
to the superior operation center in the electronic recognition system, such a center (in this 
paper, called the information-fusion center (IFC)) can perform the fusion function of the 
identification information. The fusion of identification information ensures greater stability 
of this information, i.e., resistance to accidental changes in sensor decisions. 

An ESM sensor is a passive sensor that captures incoming electromagnetic signals 
generated firstly by radar emitters mounted on air or maritime platforms. This sensor 
recognizes radar signals determining the values of their distinctive features. In this paper, 
we do not handle methods of radar signals recognition in detail. We do, however, use the 
information about these methods to identify platforms generating the signals according 
to STANAG 1241—NATO Standardization Agreement and DSmT. As previously stated, 
we are interested in three basic values of identification: friend, hostile, and neutral, as 
well as two secondary values: suspicious and assumed friendly. In addition, we assume 
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that in some situations, it is not possible to determine the identity of the emitter-carrier 
platform. To clarify this issue, we briefly describe the method used to determine the 
identification of the emitter-carrier platform that generated the captured signal. The sensor- 
recognition system is equipped with a database that can be divided into three components: 
a platform database, an emitter list, and a geopolitical list [21]. The platform database 
(PDB) contains information about platforms that can be met in the area of interest, along 
with their equipment with emitters. The emitter name list (ENL) includes all emitters 
corresponding to each platform of the PDB and contains the values of the signal distinctive 
features for each emitter. The values of distinctive features are the basis for the procedure 
of recognizing a captured signal. The geopolitical list (GPL) provides the allegiance of 
various countries and platforms and yields the identification of them in accordance with 
STANAG 1241. 

The algorithm of signal recognition is realized in two stages: 

1. Verification at the level of signal quality features. The second stage is executed after 
a positive assessment of the conformity of quality features; 

2. The signal-recognition procedure determines the distances between the distinctive 
features of the recognized signal and the distinctive features of all pattern signals stored 
within the emitter list. 

Let us introduce the following notation: 


xs : vector of distinctive features of the recognized signal; 

x;: vector of distinctive features of the i-th pattern signal (i: the number of the pattern 
signal, i € {1,...,M}); 

d,; = d(xs,x;) : the distance between the distinctive features vector of the recog- 
nized signal and the distinctive features vector of the i-th pattern signal; the dis- 
tance d,; is the Mahalanobis distance, taking into account the correlations of the 
distinctive features. 


The signal-recognition classifier compares the distance d(x,,x;) with the acceptable 
positive distance of the classification 5. The distance 6 is the limit that we interpret as a 
boundary of emitter pattern recognition. We divide the set of pattern signals into two sub- 
sets: the patterns satisfying the positive classification condition in relation to the recognized 
signal s (D{ ) and the patterns that do not satisfy the positive classification condition (D; ). 
The formal definition is as follows: 


Dy = {i€ {1,...,M}|dg; < 5} (73) 


Dy = {i€ {1,...,M}|ds; > 5} (74) 


In this paper, we propose the following method of determining the basic belief assign- 
ment on a set of pattern signals, which is related to the distance between a signal and a 
pattern in the distinctive features space: 


ms(i) = e406) (75) 


As can be seen from Formula (75), if d(xs,x;) = 0, then ms;(i) = 1, whereas if d(xs,x;) > 
0, then 0 < m;(i) < 1. The above measure is not normalized; hence, we normalize it 
as follows: 


(76) 


The sum of the measures assigned to all the emitters, the distinctive features of 
which lie outside the limit 5, are treated as a measure assigned to the base hypothesis 
“unknown” (U): 


ffis(U) = J ms(i) (77) 


i€D; 
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Another way of recognizing emitters based on their signals is presented in [26]. For 
this, the authors use a convolutional neural network with a softmax layer. 

To determine the belief measure of other base hypotheses (H, F, N) and secondary 
hypotheses (AF and S), we introduce formal definitions of sets contained in the sensor 
database and used for the recognition of captured signals. As mentioned above, the set 
of all the necessary data for platform identification can be divided into three sets: PDB, a 
platform database; ENL, an emitter name list; and GPL, a geopolitical list: 

PDB: the platform database contains information about all platforms observed in 
the area of interest, including information on all emitters mounted on each platform; we 
assume that one platform can have many emitters and the same type of emitters can be 
installed on many platforms; the PDB also contains information on the national affiliation 
of each platform; 

ENL: the emitter name list is a set of information about all recognized emitters in 
the area of interest; this set contains the mean values of the distinctive features of emitter 
signals (the so-called signal patterns) and their standard deviations; 

GPL: the geopolitical list contains base values of identification attributes (H, F, N) 
assigned to the various countries. 

We also introduce additional notations used in this paper: 

e PDBL: the list of platform numbers that are stored in the PDB; 
e PL(i): the set of numbers of platforms that have an emitter with number “7”; 


“an 


IPL(): the base identification attribute of the platform with number “j” determined on 
the basis of the information contained in the PDB and ENL (IPL(j) € { F, H, N}). 

The set of signal patterns satisfying the positive classification condition in relation to 
the recognized signal s, denoted as D,, can be divided into disjunctive subsets according 
to the values of the carrier platform identification features: 


DP =p up up up up, (78) 


D?*n Dt! =@, k #1, k,l € {F,H,N,AE,S} (79) 


Each subset of the set D7 for the base identification is defined as follows: 


Dj? = {i€ Di |Vj € PL(i) IPL(j) = F}, (80) 
Dj" = {ie Df |Vj € PL(i) IPL(j) = H} (81) 
D{N = {i € Df |Vj € PL(i) IPL(j) = N} (82) 


In a similar way, subsets of the set D for the secondary identification (AF, S) can be 
defined as follows: 


Di“" = {ie DS |Aj € PL(i) IPL(j) =F A 3j€PL(i) IPL(Gj)=N} (83) 


DiS = {ie Df :4j € PL(i) IPL(j) =H A 4j € PL(i) IPL(j)=N} (84) 


It can be noticed that we assume in this paper that no emitter type can be installed 
simultaneously on platforms with identifications F and H: 


{ieD{|4j¢PL(i) IPL(j)=F A 3j¢PL(i) IPL(j)=H}=9 (85) 


Introducing the definition of subsets of the set determines the belief masses for all 
identification features: 


ms(F)= J)" ma(t), mH) = > ms(t),ms(N) = S| m(t), (86) 


ieD?F icD}" icD{N 
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itis(AF) = Yo itis(i),iis(S) = YO fis(i) (87) 


ieD{AF ieD{$ 


It should be emphasized that the method presented here is different than that presented 
in [25,27]. These papers assume that ESM sensors can only generate basic declarations with 
attribute values FRIEND, HOSTILE, and NEUTRAL, but in this paper, we assume that 
ESM sensors can generate declarations from an extended set of attribute values (including 
ASSUME FRIEND, SUSPECT, and UNKNOWN). 


6. Numerical Experiments of Fusion of Identification Information from ESM Sensors 
6.1. General Research Scheme of Fusion of Identification Information from ESM Sensors 


Figure 4 shows a general scheme of simulation experiments, which indicates the places 
of description of individual models. 


Simulation Identification 

scenario information 
generator for ESM sensor fusion 

electronic reports carried out in 


situation tracking accordance with 


by ESM sensors one of the PCR 
described in rules described in 


Section 5 Section 3 


Figure 4. The general diagram of simulation experiments of fusion of identification information from 
ESM sensors. 


6.2. Simulation Scenarios 


Paper [25] presents a typical simulation scenario for testing identification information 
fusion. The authors formulated several requirements that should be met by such a scenario. 
It should meet the following requirements: 


(1) adequately represent the known ground truth of the emitter identification; 

(2) include sufficient numbers of incorrect associations to be realistic and to test the 
robustness of the rules in temporary incorrect sensor decisions; 

(3) provide only partial knowledge about the ESM sensor declarations and thus contain uncertainty; 

(4) be able to show stability in the case of countermeasures; 

(5) be able to switch identification when the ground truth changes. 


The authors of [25] propose the following parameters of the scenario: 


(1) ground truth of identification is FRIEND (F) for the first 50 iterations of the scenario 
and HOSTILE (H) for the last 50 iterations; 

(2) the percentage of correct associations is 80% of all iterations, and the percentage of 
incorrect associations caused by countermeasures is 20% of all iterations in randomly 
selected moments of time; 

(3) ESM sensor declarations have a mass of 0.7 for the most credible identification and 
0.3 for the identification of UNKNOWN (U). 


Assumption (5) is not considered in this paper, assuming that the real object does not 
change its real identity while performing the mission. Therefore, assumption (1) regarding 
the scenario parameters becomes obsolete. 

The following assumptions concerning the parameters of the scenario have been made 
in this paper: 

(1) the real value of identification is constant in each scenario and is equal to FRIEND (F) 

in scenarios 1, 2, and 5 and HOSTILE (H) in scenarios 3, 4, and 6; 
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(2) the above declarations are transmitted by sensor number 1 with the real identification 
mass equal to 0.7 and the mass of complementary identification (UNKNOWN) equal 
to 0.3; 

(3) the second sensor transmits its declarations in accordance with Tables 2 and 3 for 
scenarios 1 and 2 and in accordance with Tables 4 and 5 for scenarios 3 and 4. 


Table 2. Belief mass values for the second sensor for scenarios 1 and 5. 


Type of Identification F N H AF S Uu 
Correct identification (80% of events) 0.6 0.1 0 0.2 0 0.1 
Incorrect identification (20% of events) 0 0.1 0.6 0 0.2 0.1 


Table 3. Belief mass values for the second sensor for scenario 2. 


Type of Identification F N H AF S Uu 
Correct identification (80% of events) 0.7 0.1 0 0.1 0 0.1 
Incorrect identification (20% of events) 0 0.1 0.7 0 0.7 0.1 


Table 4. Belief mass values for the second sensor for scenario 3. 


Type of Identification F N H AF S Uu 
Correct identification (80% of events) 0 0.1 0.6 0 0.2 0.1 
Incorrect identification (20% of events) 0.6 0.1 0 0.2 0 0.1 


Table 5. Belief mass values for the second sensor for scenarios 4 and 6. 


Type of Identification F N H AF S Uu 
Correct identification (80% of events) 0 0.1 0.7 0 0.1 0.1 
Incorrect identification (20% of events) 0.7 0.1 0 0.1 0 0.1 


It should be noted that scenario 2 differs from scenario 1 by a greater belief mass 
assigned to an incorrect identification of the recognized emitter. Scenarios 3 and 4 are 
similarly different. 

Scenarios 1-6 for sensor 1 are presented in Figures 5 and 6. All deterministic scenarios 
for sensor 2 are presented in Figures 7-10. 


—6— Unknown 
0.9 —t— Friend 
Neutral 
0.8 —S— Hostile 
_ Assumed Friend 
O07 FER HHH] —=— Suspect HH 


mass value m(.) 
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> o 


0 10 20 30 40 50 60 70 80 90 100 
time step 


Figure 5. The course of scenarios number 1, 2, and 5 for sensor 1. 
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Figure 6. The course of scenarios number 3, 4, and 6 for sensor 1. 
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Figure 7. The course of scenario number 1 for sensor 2. 
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Figure 8. The course of scenario number 2 for sensor 2. 


This section is divided by subheadings. It provides a concise and precise description 
of the experimental results, their interpretation, and the experimental conclusions that can 
be drawn. 
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Figure 10. The course of scenario number 4 for sensor 2. 


The paper also uses the Monte Carlo method for generating the scenario for sensor 2. 
Moments in which incorrect identifications occurred are generated by the pseudorandom 
integer number generator from the range [0, 100]. Examples of scenarios are shown in 
Figures 11 and 12. 
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Figure 11. The course of Monte Carlo scenario number 5 for sensor 2. 
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Figure 12. The course of Monte Carlo scenario number 6 for sensor 2. 


6.3. Calculation Results for Deterministic Scenarios 


In all figures presenting the values of the resulting belief mass, the decision threshold 
is marked with a horizontal line. An identification whose belief mass at a given moment is 
above the decision threshold is the so-called hard decision. 


6.3.1. Dempster’s Rule 


Dempster’s rule is not resistant to a situation where the degree of conflict kp; = 1. This 
means the total conflict between the mass vector sent by the sensor and the mass vector 
of the information-fusion center, which occurs when each non-zero belief mass value sent 
by the sensor corresponds to the zero belief mass value of the vector determined by the 
information-fusion center and vice versa. 

The simulation results of identification information fusion using Dempster’s rule are 
presented for deterministic scenarios 1 and 3 in Figures 13 and 14, respectively. When 
the degree of conflict kp; = 1, according to Equation (8), it is impossible to perform sensor 
information fusion, i.e., it is impossible to determine the resulting BBA. 


‘f 
0.9 


0.8 

0.7 | 
= —8— Unknown 

0.6 
5 —t— Friend 
= 05 ——— Neutral 
> —S— Hostile 
8 ~~ Assumed Friend} | 
Py 0.4 —-A— Suspect 

0.3 

0.2 

0.1 

0 n n is 1 1 ni n 4 
0 20 40 60 80 100 120 140 160 180 200 
time step 


Figure 13. The values of the resulting belief mass for scenario 1 and Dempster’s rule. 


6.3.2. The PCR1 Rule 


The simulation results of identification information fusion using the PCR1 rule for 
deterministic scenarios 1, 2,3, and 4 are presented in Figures 15-18, respectively. 
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Figure 14. The values of the resulting belief mass for scenario 3 and Dempster’s rule. 
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Figure 15. The values of the resulting belief mass for scenario 1 and the PCRI rule. 
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Figure 16. The values of the resulting belief mass for scenario 2 and the PCRI1 rule. 
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Figure 17. The values of the resulting belief mass for scenario 3 and the PCRI1 rule. 
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Figure 18. The values of the resulting belief mass for scenario 4 and the PCR1 rule. 


6.3.3. The PCR3 Rule 


The simulation results of identification information fusion using the PCR3 rule for 
deterministic scenarios 1, 2, 3, and 4 are presented in Figures 19-22, respectively. 
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Figure 19. The values of the resulting belief mass for scenario 1 and the PCR3 rule. 
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Figure 20. The values of the resulting belief mass for scenario 2 and the PCR3 rule. 
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Figure 21. The values of the resulting belief mass for scenario 3 and the PCR3 rule. 
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Figure 22. The values of the resulting belief mass for scenario 4 and the PCR3 rule. 
6.3.4. The PCR4 Rule 


The simulation results of identification information fusion using the PCR4 rule for 
deterministic scenarios 1, 2, 3, and 4 are presented in Figures 23-26, respectively. 
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Figure 23. The values of the resulting belief mass for scenario 1 and the PCR4 rule. 
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Figure 24. The values of the resulting belief mass for scenario 2 and the PCR4 rule. 
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Figure 25. The values of the resulting belief mass for scenario 3 and the PCR4 rule. 


Remote Sens. 2023, 15, 3977 26 of 45 


Assumed Friend 
—4— Suspect 


mass value m(.) 
° 
a 


time step 


Figure 26. The values of the resulting belief mass for scenario 4 and the PCR4 rule. 


6.3.5. The PCR5 Rule for 2 BBAs 


The simulation results of identification information fusion using the PCRS rule for two 
BBAs for deterministic scenarios 1, 2, 3, and 4 are presented in Figures 27-30, respectively. 
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Figure 27. The values of the resulting belief mass for scenario 1 and the PCR5 rule for 2 BBAs. 
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Figure 28. The values of the resulting belief mass for scenario 2 and the PCR5 rule for 2 BBAs. 
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Figure 29. The values of the resulting belief mass for scenario 3 and the PCR5 rule for 2 BBAs. 
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Figure 30. The values of the resulting belief mass for scenario 4 and the PCR5 rule for 2 BBAs. 


6.3.6. The PCR5 Rule for 3 BBAs 


The simulation results of identification information fusion using the PCRS rule for three BBAs 
for deterministic scenarios 1, 2,3, and 4 are presented in Figures 31-34, respectively. 
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Figure 31. The values of the resulting belief mass for scenario 1 and the PCR5 rule for 3 BBAs. 
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Figure 32. The values of the resulting belief mass for scenario 2 and the PCR5 rule for 3 BBAs. 
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Figure 33. The values of the resulting belief mass for scenario 3 and the PCR5 rule for 3 BBAs. 
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Figure 34. The values of the resulting belief mass for scenario 4 and the PCR5 rule for 3 BBAs. 


6.3.7. The PCR6 Rule for 3 BBAs 


The simulation results of identification information fusion using the PCR6 rule for three 
BBAs for deterministic scenarios 1, 2, 3, and 4 are presented in Figures 35-38, respectively. 
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Figure 35. The values of the resulting belief mass for scenario 1 and the PCR6 rule for 3 BBAs. 
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Figure 36. The values of the resulting belief mass for scenario 2 and the PCR6 rule for 3 BBAs. 
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Figure 37. The values of the resulting belief mass for scenario 3 and the PCR6 rule for 3 BBAs. 
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Figure 38. The values of the resulting belief mass for scenario 4 and the PCR6 rule for 3 BBAs. 


The presented results (Figures 13-38) yield the conclusion that the applied methods of 
managing conflicts in information fusion enables correct conclusions to be drawn about the 
real identification of the recognized object. 

The application of the decision threshold for the belief mass at the level my = 0.37 for 
the PCR1 rule (Figures 15 and 17) and my, = 0.45 for PCR3, PCR4 (Figures 19, 21, 23 and 25), 
and PCR5 for two BBAs (Figures 27 and 29) for scenarios 1 and 3 allows for a proper 
evaluation of the identification of the recognized object: scenario 1, FRIEND; scenario 3, 
HOSTILE. For scenarios 2 and 4, the optimal thresholds are mg = 0.4 for the PCR1 rule 
(Figures 16 and 18) and my, = 0.48 for the PCR3, PCR4 (Figures 20, 22, 24 and 26), and 
PCR5 rules for two BBAs (Figures 28 and 30). When assessing the interval between the 
minimum resultant mass for correct identification and the maximum resultant mass for 
misidentification, the worst results are reached by the PCR1 rule, and the rules of PCR3, 
PCR4, and PCR5 behave similarly and are better than rule PCR1. 

The research carried out for the deterministic scenarios shows that the PCR5 rule for 
three BBAs and the PCR6 rule for three BBAs behave very similarly (Figures 31, 33, 35 and 37 
for scenarios 1 and 3, mg = 0.37 and Figures 32, 34, 36 and 37 for scenarios 2 and 4, my = 0.39). 
They restore the correct identification after the occurrence of temporary misidentification 
much faster than the rules PCR1-PCR5 for two BBAs. 


6.4. Calculation Results for the Monte Carlo Scenarios 
6.4.1. Dempster’s Rule 


In the Monte Carlo scenario, Dempster’s rule behaves similarly to a deterministic 
scenario. It is not resistant to a situation where the degree of conflict kp; = 1. This means 
that the total conflict between the mass vector sent by the sensor and the mass vector of 
the information-fusion center, which occurs when each non-zero belief mass value sent 
by the sensor corresponds to the zero belief mass value of the vector determined by the 
information-fusion center and vice versa. 

The simulation results of identification information fusion using Dempster’s rule are 
presented for Monte Carlo scenarios 5 and 6 in Figures 39 and 40, respectively. 
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Figure 39. The values of the resulting belief mass for Monte Carlo scenario 5 and Dempster’s rule. 
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Figure 40. The values of the resulting belief mass for Monte Carlo scenario 6 and Dempster’s rule. 


6.4.2. The PCR1 Rule 


The simulation results of identification information fusion using the PCR1 rule for 
Monte Carlo scenarios 5 and 6 are presented in Figures 41 and 42, respectively. The 
application of the decision threshold for the belief mass at the level mg, = 0.34 for the 
PCR1 rule (Figures 41 and 42) for scenarios 5 and 6 allows for a proper evaluation of the 
identification of the recognized object. There are only three time points when the rule 
misidentifies due to increased misidentification intensity. 
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Figure 41. The values of the resulting belief mass for Monte Carlo scenario 5 and the PCR1 rule. 
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Figure 42. The values of the resulting belief mass for Monte Carlo scenario 6 and the PCR1 rule. 


6.4.3. The PCR3 Rule 


The simulation results of identification information fusion using the PCR3 rule for 
Monte Carlo scenarios 5 and 6 are presented in Figures 43 and 44, respectively. The 
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application of the decision threshold for the belief mass at the level mg, = 0.42 for the 
PCR3 rule (Figures 43 and 44) for scenarios 5 and 6 allows for a proper evaluation of the 
identification of the recognized object. There is only one time point for scenario 6 where 
the rule misidentifies due to increased misidentification intensity. 
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Figure 43. The values of the resulting belief mass for Monte Carlo scenario 5 and the PCR3 rule. 
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Figure 44. The values of the resulting belief mass for Monte Carlo scenario 6 and the PCR3 rule. 


6.4.4. The PCR4 Rule 


The simulation results of identification information fusion using the PCR4 rule for 
Monte Carlo scenarios 5 and 6 are presented in Figures 45 and 46, respectively. The 
application of the decision threshold for the belief mass at the level mg, = 0.47 for the 
PCR4 rule (Figures 45 and 46) for scenarios 5 and 6 allows for a proper evaluation of the 
identification of the recognized object. 


NT _—3— unknown 
—— Friend 
Neutral 
—S— Hostile 
~ Assumed Friend 
—+— Suspect 


mass value m(.) 


0 20 40 60 80 100 120 140 160 180 200 
time step 


Figure 45. The values of the resulting belief mass for Monte Carlo scenario 5 and the PCR4 rule. 
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Figure 46. The values of the resulting belief mass for Monte Carlo scenario 6 and the PCR4 rule. 
6.4.5. The PCR5 Rule for 2 BBAs 


The simulation results of identification information fusion using the PCR5 rule for two 
BBAs for Monte Carlo scenarios 5 and 6 are presented in Figures 47 and 48, respectively. 
The application of the decision threshold for the belief mass at the level 1, = 0.47 for the 
PCRS rule for two BBAs (Figures 47 and 48) for scenarios 5 and 6 allows for a proper 
evaluation of the identification of the recognized object. There is only one time point for 
scenario 6 where the rule misidentifies due to increased misidentification intensity. 
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Figure 47. The values of the resulting belief mass for Monte Carlo scenario 5 and the PCR5 rule for 
2 BBAs. 
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Figure 48. The values of the resulting belief mass for Monte Carlo scenario 6 and the PCR5 rule for 
2 BBAs. 
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6.4.6. The PCR5 Rule for 3 BBAs 


The simulation results of identification information fusion using the PCR5 rule for 
three BBAs for Monte Carlo scenarios 5 and 6 are presented in Figures 49 and 50, respec- 
tively. The application of the decision threshold for the belief mass at the level my = 0.39 for 
the PCR5 rule for three BBAs (Figures 49 and 50) for scenarios 5 and 6 allows for a proper 
evaluation of the identification of the recognized object. 
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Figure 49. The values of the resulting belief mass for Monte Carlo scenario 5 and the PCR5 rule for 
3 BBAs. 
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Figure 50. The values of the resulting belief mass for Monte Carlo scenario 6 and the PCR5 rule for 
3 BBAs. 


6.4.7. The PCR6 Rule for 3 BBAs 


The simulation results of identification information fusion using the PCR6 rule for 
three BBAs for Monte Carlo scenarios 5 and 6 are presented in Figures 51 and 52. The 
application of the decision threshold for the belief mass at the level my = 0.39 for the PCR6 
rule for three BBAs (Figures 51 and 52) for scenarios 5 and 6 allows for a proper evaluation 
of the identification of the recognized object. 
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Figure 51. The values of the resulting belief mass for Monte Carlo scenario 5 and the PCR6 rule for 
3 BBAs. 
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Figure 52. The values of the resulting belief mass for Monte Carlo scenario 6 and the PCR6 rule for 
3 BBAs. 


The presented results show that due to the high intensity of sending reports with 
incorrect identifications in the middle part of the scenarios, the information-fusion rules 
(apart from the PCR4, PCR5, and PCR6 rules) determine the maximum resulting mass for 
incorrect identification. The PCR5 for three BBAs and PCR6 for three BBAs rules are the 
fastest to restore the correct identification after receiving several incorrect reports. 


7. Numerical Experiments of Fusion of Identification Information from ESM Sensors 
and Radars 

7.1. General Research Scheme of Fusion of Identification Information from Radars and 

ESM Sensors 


Figure 53 shows a general scheme of simulation experiments, which indicates the 
places of description of individual models. 
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Figure 53. The general diagram of simulation experiments of fusion of identification information 
from radars and ESM sensors. 


7.2. Numerical Experiments Scenarios 

We assume that we will combine attribute information from two sensors: a combined 
primary and secondary surveillance radar and ESM sensor. These sensors work asyn- 
chronously. Upon receipt of the sensor’s declaration in the form of a vector of masses, we 


Remote Sens. 2023, 15, 3977 


35 of 45 


fuse this vector with the vector of the actual values of the declaration masses for the fuser’s 
frame of discernment. The frequency of transmission of the sensor declarations depends on 
the rules of the data exchange network and on the technical characteristics of the sensors. 
Various combination methods are presented in [3,4]. This paper used two of the methods 
of proportional redistribution conflict (PRC5 and PCR6 [6,22]). Information fusion has 
been simulated for two processing schemes (Figures 2 and 3). The numerical model of 
combined primary and secondary surveillance radars was taken from [17,23]. It allows for 
the determination of the probability Py during the simulation of the object’s movement, i.e., 
the change in the object’s position relative to the radar. Detailed rules for determining BBAs’ 
vectors, assuming the knowledge of probabilities Py and Pfr, are presented in Section 4. 
Numerical experiments have been performed for the following data: 


e for combined primary and secondary surveillance radars sensor: 
Prg= 10~°, Ringy= 100 [km], Pj= 0.7, of =2 |m?|, Prep = 0.962 
and the following table of masses (compare Table 6): 


Table 6. Transformation of the base belief assignment mass into the secondary belief assignment 
mass for combined primary and secondary surveillance radar. 


(Scenario Nr, Base Identification) > (1,FR) (2,Npg) (3,Hp) 
Fa 0.8 0 0 
Ns 0 0.5 0 
ae 0 0 0.7 
AF 0.2 0.3 0 
S 0 0.2 0.3 


The flight path of the air object was 30 km away from the sensor (in the horizontal 
plane), the flight altitude was 1 km, and the radar cross-section was 1 m2. 

The following assumptions concerning the parameters of the scenario for the ESM 
sensor were made in this paper: 

(1) The real value of identification is constant in each scenario and is equal to FRIEND 
(F) in the first scenario and HOSTILE (H) in the second scenario; 

(2) The above declarations are transmitted by sensor number 1 with the real identifica- 
tion mass equal to 0.7 and the mass of complementary identification (UNKNOWN) equal 
to 0.3; 

(3) The second sensor shall transmit its declarations in accordance with Tables 1 and 2 
for scenarios 1 and 2, respectively, and with Tables 3 and 4 for scenarios 3 and 4, respectively. 


Tables 7-12 present the mass values for all possible declarations for the six scenarios 
for the ESM sensor. 


Table 7. Belief mass values for the second sensor (ESM) for scenario 1. 


Type of Identification F N H AF S Uu 
Correct identification (80% of events) 0.6 0.1 0 0.2 0 0.1 
Incorrect identification (20% of events) 0 0.1 0.6 0 0.2 0.1 


Table 8. Belief mass values for the second sensor (ESM) for scenario 2. 


Type of Identification F N H AF S U 


Correct identification (80% of events) 0 0.5 0.3 0.2 0 
Incorrect identification (20% of events) 0 0.4 0.2 0 0.3 0.1 
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Table 9. Belief mass values for the second sensor (ESM) for scenario 3. 

Type of Identification F N H AF S U 

Correct identification (80% of events) 0 0.1 0.7 0 0.1 0.1 

Incorrect identification (20% of events) 0 0.1 0.6 0 0.2 0.1 
Table 10. Belief mass values for the second sensor (ESM) for scenario 4. 

Type of Identification F N H AF S U 

Correct identification (80% of events) 0.1 0.7 0.1 0 0 0.1 

Incorrect identification (20% of events) 0 0.1 0.6 0 0.2 0.1 
Table 11. Belief mass values for the second sensor (ESM) for scenario 5. 

Type of Identification F N H AF S U 

Correct identification (80% of events) 0.6 0.1 0 0.2 0 0.1 

Incorrect identification (20% of events) 0 0.1 0.6 0 0.2 0.1 
Table 12. Belief mass values for the second sensor (ESM) for scenario 6. 

Type of Identification F N H AF S U 

Correct identification (80% of events) 0.1 0.7 0.1 0 0 0.1 

Incorrect identification (20% of events) 0.6 0.1 0 0.2 0 0.1 


Scenarios 1-6 for sensor 1 are presented in Figures 54-56. All deterministic scenarios 


for sensor 2 are presented in Figures 57-62. 
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Figure 54. The course of scenarios number 1 and 4 for sensor 1. 
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Figure 55. The course of scenarios number 2 and 5 for sensor 1. 
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Figure 56. The course of scenarios number 3 and 6 for sensor 1. 
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Figure 57. The course of scenario number 1 for sensor 2. 
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Figure 58. The course of scenario number 2 for sensor 2. 
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Figure 59. The course of scenario number 3 for sensor 2. 
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Figure 60. The course of scenario number 4 for sensor 2. 
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Figure 61. The course of scenario number 5 for sensor 2. 
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Figure 62. The course of scenario number 6 for sensor 2. 
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Scenarios 1-3 assume relatively small changes in the mass of all declarations. Scenarios 1-3 
assume significant changes in the credibility mass of all declarations (small errors). Scenarios 4-6 


assume significant changes in the mass of all declarations (large errors). 


7.3. Calculation Results for Four Proportional Conflict Redistribution Rules 


7.3.1. The PCR5 Rule for 2 BBAs 


The simulation results of identification information fusion using the PCR5 rule for 


two BBAs for deterministic scenarios 1-6 are presented in Figures 63-68. 
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Figure 63. The values of the resulting belief mass for scenario 1 and the PCR5 rule for 2 BBAs. 
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Figure 64. The values of the resulting belief mass for scenario 2 and the PCR5 rule for 2 BBAs. 
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Figure 65. The values of the resulting belief mass for scenario 3 and the PCR5 rule for 2 BBAs. 
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Figure 66. The values of the resulting belief mass for scenario 4 and the PCR5 rule for 2 BBAs. 
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Figure 67. The values of the resulting belief mass for scenario 5 and the PCR5 rule for 2 BBAs. 
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Figure 68. The values of the resulting belief mass for scenario 6 and the PCR5 rule for 2 BBAs. 
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Figure 69. The values of the resulting belief mass for scenario 1 and the PCR5 rule for 3 BBAs. 
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Figure 70. The values of the resulting belief mass for scenario 2 and the PCR5 rule for 3 BBAs. 
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Figure 71. The values of the resulting belief mass for scenario 3 and the PCR5 rule for 3 BBAs. 
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Figure 72. The values of the resulting belief mass for scenario 4 and the PCR5 rule for 3 BBAs. 
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Figure 73. The values of the resulting belief mass for scenario 5 and the PCR5 rule for 3 BBAs. 
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Figure 74. The values of the resulting belief mass for scenario 6 and the PCR5 rule for 3 BBAs. 
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The simulation results of identification information fusion using the PCR6 rule for 
three BBAs for deterministic scenarios 1-6 are presented in Figures 75-80. 


Figure 75. The values of the resulting belief mass for scenario 1 and the PCR6 rule for 3 BBAs. 
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Figure 76. The values of the resulting belief mass for scenario 2 and the PCR6 rule for 3 BBAs. 
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Figure 77. The values of the resulting belief mass for scenario 3 and the PCR6 rule for 3 BBAs. 
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Figure 78. The values of the resulting belief mass for scenario 4 and the PCR6 rule for 3 BBAs. 
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Figure 79. The values of the resulting belief mass for scenario 5 and the PCR6 rule for 3 BBAs. 
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Figure 80. The values of the resulting belief mass for scenario 6 and the PCR6 rule for 3 BBAs. 
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For the PCR6 rule for three BBAs, the application of the decision thresholds at the 
belief mass level m, = 0.45 for scenario 4 allows for a proper evaluation of the identification 
of the recognized object. 

For the PCR6 rule for three BBAs, the application of the decision thresholds at the 
belief mass level m, = 0.35 for scenario 5 and m, = 0.4 for scenario 6 allows for a proper 
evaluation of the identification of the recognized object for most time moments. 

Comparing Figures 66-68 with Figures 72-74 and 78-80, conclusion can be drawn 
that the PCR5 for three BBAs and PCR6 for three BBAs rules provide more stable results 
of combined belief masses (smaller amplitude of changes). Due to the large dispersion 
of belief mass changes for scenarios 5 and 6, it is not possible to correctly evaluate the 
identification of the recognized object for all time moments. 

The presented results (Figures 61-78) allow a conclusion to be drawn that the applied 
methods of removing conflicts in information fusion enables the correct conclusions to be 
drawn about the real identification of the recognized object. 


8. Conclusions 


The proposed basic belief assignment model for ESM sensors and radars can be used 
to build identification information-fusion systems. Models conformable to STANAG 1241 
have primary practical significance. 

Due to the assumption of conflicts between the ESM sensor declarations in this work, 
Dezert-Smarandache theory is used to determine the basic belief assignment of decla- 
rations as a product of the process of fusion of identification information sent by these 
sensors. Supplementing standard reports on the detected signals with random identifi- 
cation declarations allows the use of methods of identification information fusion in the 
information-fusion center. The test results confirm the full usefulness of conflict redistribu- 
tion rules in reports from ESM sensors developed as a part of Dezert-Smarandache theory, 
with the best results presented for the PCR5 and PCR6 rule. 

The basic belief assignment model for ESM sensors and for combined primary and 
secondary radars [17] can be applied to build models of different identification data-fusion 
systems. All the models compatible with STANAG 1241 have primary practical significance 
as it contains definitions corresponding to intersections of basic identification declarations. 
Therefore, the paper uses Dezert-Smarandache theory for the calculation of the basic 
belief assignment. 

The conducted research showed that the best results were obtained for the PCR6 rule 
when reports from three sources (from two sensors and the fusion-system database) were 
processed simultaneously. This corresponds to the synchronous processing of reports and 
involves delayed processing of a report from one of the sources. The research confirmed a 
slight advantage of the PCR6 rule over the PCR5 rule. This was mainly the case when the 
sensors sent information with a high degree of conflict. 
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